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1. Introduction 

Polymers and polymer solutions belong to the most intensively studied objects in 
condensed matter physics [H 121 El- The adsorption of polymers on surfaces and interfaces 
is of special importance jlj. Adsorption of free polymers in solution on the container 
wall or other boundaries occurs in the presence of attractive interactions between the 
surface and the chain monomers. Examples of such interactions include Coulomb and 
van der Waals forces, and more specific molecular interactions. Coulomb forces are 
screened by counter-ions in the solution and can be tuned to some extent by adding 
salt to the solvent. Thus on changing the properties of the solvent an individual 
polymer chain can undergo a reversible transition from a freely floating, desorbed 
state to an adsorbed state in which the chain monomers are close to the surface on 
average. The adsorption of polymers on flat surfaces has been studied theoretically and 
experimentally, and is by now well understood [TJ HJ El Ej. Due to the importance to 
colloidal dispersions the interaction of polymers with spherical and rodlike particles has 
been studied as well [11 El HE EDI HI] • The adsorption of flexible polymers on rodlike 
particles is relevant, for example, in gels [Oj, and for the binding of flexible polymers 
to relatively stiff biomolecules such as DNA ^2] • Another class of polymer adsorption 
transitions involve two flexible self-avoiding but mutually attracting polymers which can 
form a bound, double-stranded, so-called "zipped" state. A prominent example of this 
kind of transition is the denaturation transition of double-stranded DNA [T31 ITlj which 
recently attracted considerable attention regarding its statistical-mechanical properties 
and thermodynamic order [To^ \W\ ITT| HHl E3 120] (for a recent review on the biophysics 
related to DNA topology, see [21] )• The DNA denaturation transition is usually 
modelled in such a way that monomer % of one strand can only interact with monomer i of 
the other strand, reflecting the key-lock principle of natural, inhomogeneous DNA with 
its specific sequence of base pairs. Two-chain systems in which any monomer of one chain 
can interact with any monomer of the other chain include diblock copolymers, which 
consist of linear chains of N monomers of type A followed by N' monomers of type B, 
with different AA, BB, and AB interactions; systems of this kind have been intensively 
studied as well [22 E3 I2H 123 I2E1 121] • Recently it was found that self-avoiding mutually 
attracting diblock copolymers can adopt a zipped state in which the two components 
form a bound, double-stranded structure which however remains swollen and does not 
assume compact configurations. The zipped state is located between a swollen, unbound 
high-temperature state and a compact low-temperature state, separated by second-order 
and first-order phase transitions, respectively j2Hl EH1 ED] • 

In this work we focus on the adsorption transition of a long flexible self-avoiding 
polymer chain onto a rigid thin elongated rod, as shown in figure H We assume that 
the rod is endowed with a short-ranged surface potential, or adsorption energy, £ , for 
the chain monomers; the more positive £, the more attractive the interaction. Thus, 
on increasing S from a low value, at some threshold value £* the chain is expected to 
undergo a transition from an unbound, free state to a bound state in which the overall 
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Figure 1. A long flexible self- avoiding polymer chain (green) adsorbs onto a thin 
elongated rod (black) if the adsorption energy per chain monomer £ is larger than the 
threshold value £* . 



gain in binding energy compensates the loss of configurational entropy. An interesting 
feature of this transition is the fact that it represents a true phase transition in the 
thermodynamic sense; the adsorbed state forming an elongated, aligned structure, see 
figure H In contrast, for a particle of finite size, a self-avoiding polymer chain cannot 
undergo a true adsorption phase transition due to steric constraints. The focus of 
the present work is on the former case. Since the polymer adsorption transition is 
characterized by large fluctuations we expect scaling and universal behaviour. We 
thus use the renormalization group to obtain general scaling results for the chain 
partition function. We also obtain estimates for relevant critical exponents by geometric 
arguments. Before we introduce and study our model in section El we recall some 
general ideas and concepts for later reference. In section 11.11 we discuss typical scaling 
arguments related to the polymer adsorption transition. Since our work strongly relies 
on field-theoretical methods, in section 11.21 and | Appendix A we sketch the mapping of 
the polymer system on the Ginzburg-Landau model. 



1.1. General scaling behaviour 

Consider the adsorption of a long flexible polymer chain onto an object S. For the time 
being, this object can be a surface, a thin rod, another flexible polymer chain, or any 
other extended manifold that allows for a thermodynamic adsorption phase transition. 
The quantity of interest is the partition function Z of the chain with one end fixed 
close to the object S and the other end moving freely. Upon adjusting the system's 
thermodynamic variables close to the adsorption transition point, only the number 
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of chain monomers N ^> 1 and the adsorption energy 8 ~ 8* remain as relevant 
parameters, where 8* is the adsorption energy at the transition point. The partition 
function Z is expected to obey the scaling 

Z(N, 8) ~ p N N^- l f{{8 - 8*)N*} (1) 

where p is the lattice-dependent connectivity constant and 7', $ are critical exponents. 
The scaling function f(x) is regular at x = since Z(N, 8) has no singularity for finite N 
and 8 ~ 8*. The exponent 7' thus characterizes the scaling of Z right at the transition 
point: Z(N,8*) ~p N N~t'-\ Note that 7' is not necessarily equal to the critical exponent 
7 for an unbounded, free chain, for which Zf ree (N) ~ p N N 1 ^ 1 (compare equation (|56|) in 
section l2~3t with 71 introduced in equation (}52|) and L ~ iV). The exponent <£> is referred 
to as the crossover-exponent. Since —8 acts as a chemical potential for monomers close 
to S, the number N s of adsorbed monomers scales as 

±\nZ{N,8) . (2) 

Equation ([!} implies three distinct scaling regimes for Ng. 

(i) 8 = 8*. Equations (0) and © yield 

iVs ~ AT* , 8 = 8* , N -> 00 . (3) 

For < $ < 1 this implies that iVs grows with iV but the fraction of adsorbed monomers, 
Ns/N, vanishes for N — > 00. For $ = 1, the behaviour iVg ~ at £^ = indicates 
that the adsorption transition is of first order. 

(ii) 8 < 8*. Equation (0) implies that the scaling behaviour of Z for iV — > cxd is 
governed by the behaviour of /(x) for x — > —00, regardless the precise value of £. In 
this case 5 is repulsive for the chain monomers and Ns stays finite even for iV — > 00. 

(iii) 8 > 8*. The chain adopts an adsorbed state and stays close to S on average. 
Thus, Ns ~ AT, which implies a finite fraction of adsorbed monomers for A^ — > 00: 

^ hm > , 5 > r . (4) 

To analyze the behaviour of it is useful to consider the grand canonical ensemble. 
The partition function in the grand canonical ensemble, X(^i,8), is related to Z(N,8) 
by a Laplace transform: 

POO 

X(p,£)= dNe-^ 'Z(N,8) , (5) 







where p, is the chemical potential conjugate to N. Equation © is valid for p > p c with 
p c = In p. One is allowed to set p = 1 for simplicity, so that p c = 0. Equation JIJ then 
implies the scaling behaviour 

X(ji, 8) ~ - /i > . (6) 

By reasoning similar to that below equation $1]) one finds that the scaling function g(y) 
is regular at y = and 7' characterizes the scaling of <Y right at the transition point: 
X(p,8*) ~ pT 1 ' . On the other hand, we know that for 8 > 8* the chain takes an 
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Figure 2. Typical phase diagrams for polymer adsorption in the a) grand canonical 
ensemble (/x, £) and b) canonical ensemble (N,£). The manifold S on which the 
polymer adsorbs can be any object which allows for an adsorbed thermodynamic state. 
For N — > co, the fraction of adsorbed monomers N$/N vanishes for £ < £ * (blue lines) 
whereas it is finite for £ > £* (red lines). The shape of the red line in a) is given by 
the function /is(£) in equation JHJ. Compare figure ^ 



adsorbed state. For the grand canonical ensemble this implies, for given £ > £*, that 
Ns — > oo for fM \ fis(£) with some fis(£) > 0. In this limit we thus expect the scaling 
behaviour X(fi,£) ~ (/x — /is 1 ) -75 where 75 is characteristic for the adsorbed state and 
different from 7'. For example, if S is another flexible polymer chain, the adsorbed state 
forms a double-stranded structure which, as a whole, behaves like an unbounded, free, 
self-avoiding chain, which implies 75 = 7 in this case pH] . Using equation (jHJ) it follows 
that the scaling function g(y) must have a singularity at some y$ > of the form 

g(y /ys)~{ys-yr is • (7) 

The relation {£ — £*)^i s ^ = ys determines the shape of the line Hs(£) as 

Hs{£) ~ {£ - ,£>£*■ (8) 

Figure El shows typical phase diagrams for polymer adsorption in the grand canonical 
ensemble (fixed fi) and canonical ensemble (fixed N). 

According to the above, for £ > £*, X(fi,£) in equation (0) has a singularity for 
H \ Hs(£) > 0, and N (> Ns) diverges in this limit. This, in turn, implies 

Z(N,£) ~ exp[/i 5 (£)iV] , £ > £* , iV -> oo . (9) 

Using equation (J2J one finds 

N s ~N-jtii 8 {£) , (10) 
and thus, using equations (jlj) and (jHJ), 

F(£)= lhn ~ ± p 8 (£) ~ (£ - £•)* , (11) 
n^oo Jy d£ 

where the exponent k is related to the crossover-exponent $ in equation Q by 
1 - $ 

« = -|-- (12) 
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In particular, for $ = 1 the fraction lim N$/N jumps from zero for £ < £* to a finite 

N-^QO 

value for £ > £*, which is then independent of £; this corresponds to a first-order 
transition (compare case (i) above). 

The scaling behaviours (JHJ) and (JT2|) have been proven rigorously for the 

polymer adsorption transition on a flat surface IE]- This system is closely related 
to the semi-infinite Ginzburg-Landau model, see jSH E21 EB] for reviews; the mapping 
of the polymer system on the Ginzburg-Landau model is discussed below. The scaling 
behaviours (J3J) and (JHJ), (JT2"J) also hold reasonably well in a recent numerical study of 
the DNA denaturation transition [TK] . 

1.2. Mapping of the polymer system on the Ginzburg-Landau model 

According to Edwards' continuous chain model we represent the configuration of a linear 
chain of length L by a curve R(s), parameterized by its arc length s, in D-dimensional 
space. The chain length L is proportional to the number of chain monomers N. In the 
presence of an external potential V(r) the partition function of the chain is given by 

-L / Jn \ 2' 

(13) 



x exp { — I d D r 



\/(r)p(r) + ^ 2 (r)]} 



with the monomer density 



L 



p(r) = / dsd u (r-K(s)) . (14) 
Jo 

rr' 

The symbol J r T>R denotes functional integration over all chain configurations with 
the chain ends fixed at r and r'. The superscript "(2)" on Z indicates that the chain 
is fixed with both ends. The coupling constant u of the p 2 (r) interaction characterizes 
the strength of the contact interaction between chain monomers: u = describes a 
Gaussian random walk whereas u > describes a self-avoiding chain. The case u < is 
related to the polymer collapse transition to a compact state in a poor solvent [Ml ESI; 
in this work we do not consider this collapse transition, and therefore exclude the case 
u < 0. 

As first noticed by de Gennes [36J, the polymer system can be mapped on the 
Ginzburg-Landau model of an n-component order parameter field $ = ($i, . . . , $ n ) 
in the limit n — > 0. It is worthwhile to note that this mapping not only works in 
perturbation theory but already on the level of the Hamiltonian in the Ginzburg-Landau 
model. The derivation, using a Gaussian transformation to linearize the p 2 (r) interaction 
in equation (|T5j) [3*7] . is left to Appendix A| The result is 



Z^(r,r';L) = C t ^ L lim (^(r) ^(r')) (15) 

where 

($i(r) $i(r')> = J Vi $i(r) ^(r') e - n{$} (16) 
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is the two-point correlation function in the Ginzburg-Landau model with the standard 
Hamiltonian 



H{$} = J al D r 
The operation 



I(V$) 2 + ^ 2 + iv(r)$ 2 + ^($ 2 ) 2 



(17) 



Ct^ L = ^~ / dte tL (18) 

acting on the correlation function in equation (|15|) is an inverse Laplace transform in 
which the integration path C in the complex t-plane is a parallel to the imaginary axis 
to the right of all singularities. Equations (J15)) - (|18|) describe the statistics of the 
polymer chain in terms of properties of the near-critical ferromagnetic n-vector model 
in the limit n — > 0. In the context of the n-vector model, the parameter t ~ T — T c 
describes the deviation of the temperature T from the critical temperature T c . The 
form of the interaction involving the potential V(r) in equation (fTTjl shows that the 
0(n)-invariant scalar <3> 2 is related to the monomer density p(r) of the polymer chain 
in equation ()13J) . Thus, translated to the polymer system, the term $ 2 • $ 2 is related 
to the contact interaction between chain monomers and t plays the role of a chemical 
potential for chain monomers in the bulk. 



2. Polymer adsorption transition onto a rigid rod 

The objective of this work is the study of the adsorption of a long flexible polymer chain 
onto a rigid thin rod. However, for the time being we model the rod by an infinitely 
elongated cylinder with small but finite radius R. The introduction of a finite cylinder 
radius R is necessary since the limit of a rod with zero radius turns out to be too singular 
for the present field-theoretical treatment; see figure El and the discussion below equation 
(I29|) . On the other hand, the fact that the adsorption transition now takes place on a 
surface, albeit a curved one, allows us to take advantage of available results for field 
theories with curved boundaries [HI] ; see section 12 .HI below. The chain of total length 
Lq ~ N, where iV is the number of chain monomers, is fixed with one end at the point 
r$ on the cylinder surface S while the other end is moving freely. The cylinder surface 
is endowed with a short-ranged surface potential Cq acting on the chain monomers, and 
it is understood that the chain monomers are excluded from the interior of the cylinder. 
The potential V(r) in equation (|T3|) is thus given by 

V(r)=c [ dS'5 D (r-r' s ) (19) 
Js 

and V(r) = oo if r is located in the interior of the cylinder. By virtue of equation (fTB^l . 
the chain partition function is given by 

Z(Lq) = [ d D r' Z^{v s , r'; L ) = C to _> Lo lim x {t ) . (20) 
Jv n ~ >0 

The integration volume V is the outer space of the cylinder bounded by the cylinder 
surface S. On the right hand side, x(to) = x{ r s]to), where the susceptibility x(r;i ) 
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rod 



curved surface 



Figure 3. The rod is represented by a cylinder with finite radius R. This allows one 
to use available results for field theories with curved boundaries, but introduces the 
additional variable R into the problem. Compare figure ^ 



in the Ginzburg-Landau model is obtained by integrating the two-point correlation 
function, i.e., 



Jv 

in a Ginzburg-Landau type field theory with Hamiltonian [2H E2 EH] 



H{<5>} 



d D r 



(V$) : 



+ ^ ($ 2 ) 2 
24 v ; 



+ 



2 



(21) 



(22) 



of an n-component order parameter field $ = ($i, . . . , $ n ). In the following we will 
understand x(r; t ) as the chain partition function in the grand canonical ensemble, 
where t is conjugate to L (compare section IT~TI with t ~ n, L ~ N, and c ~ —£)• 
The subscript "0" on Co, to, -^o, is used to distinguish these parameters from their 
renormalized counterparts that will be introduced below. If r = r,s in x(r;£o) we 
suppress the argument rs since x does not depend of the choice of by symmetry; 
compare equation (|2T)|) . 

Before we proceed with the renormalization of the model defined by equation ()22j) 
we review some cases in which results are available. To this end it is useful to consider not 
only a cylinder in D = 3 dimensions but bodies of more general shape. The "generalized 
cylinders" have an infinitely extended axis of dimension D — d and a curved surface with 
constant curvature radius R in the subspace of co-dimension d. f The axis can be the 

| The co-dimension of an object of dimension D tj in a space of dimension D is given by D — T) bj- 
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d co-dimension 




space dimension 



Figure 4. Diagram of "generalized cylinders" which behave as relevant or irrelevant 
perturbations for long, flexible, self-avoiding polymer chains. The parameter d < D 
characterizes the shape of the body and D is the space dimension. The point 
(d, D) — (2, 2) corresponds to a disc in D = 2 and the points (3, 3) to a sphere in D = 3. 
The point (2, 3) corresponds to an infinitely elongated cylinder in D = 3 (in which we 
are ultimately interested). The line D = 4 indicates the upper critical dimension 
where the polymer chains behave like Gaussian random walks. The open circles and 
the dashed line (blue) indicate points (d, D) for which the dimension D int = v~ 1 {D) — d 
of intersection points of the body with an unbounded, free, self-avoiding random walk 
vanishes. For Di nt < the intersection probability is zero. In this sense, in the shaded 
region above the dashed line the body represents an irrelevant perturbation. 

axis of an ordinary cylinder in three dimensions, for which (d, D) = (2, 3), the midplane 
of a plate (d = 1), or the centre of a sphere (d = D). The explicit form of a "generalized 
cylinder" is given by the set {r = (rj_, ry) G M, d x |r ± | < R\. Figure 0] shows some 

examples in the (d, -D)-plane; compare also references f01EJII2]- 
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2.1. Flat surface 

First, consider the line d = 1 with D arbitrary in figure El Then, of course, equation 
f!22j) corresponds to the semi-infinite n- vector model bounded by a flat surface for which 
many results are available [HU EH1 • I n particular, for the polymer case n = 0, the 
considerations of section ITTT1 can be made explicit, proving equations (jHJ) and (fT2"|) 

w 



2.2. Gaussian model 

Next, consider the line D = 4 with d arbitrary in figure EJ Since D = 4 is the upper 
critical dimension of the model in equation (|2*2*|) . this case corresponds to the Gaussian 
field theory and to Gaussian random walks, respectively. Both cases are described by 
Mo = in equation (j22J) and can be solved exactly by standard methods [12] . Thus, the 
susceptibility defined in equation (|2"Tj). corresponding to the chain partition function in 
the grand canonical ensemble, is given by ^0] 

Co P - a K a ( Py /r Q ) 



^ 1 

X(r;to) = — 
to 



1 



(23) 



/To~K a+1 (v/Tb) + Co^a (V^) 

where p = \r±\/R (so that p = 1 for r = r$), Co = Rcq and r = R 2 t . The functions K a 
and K a+ i are modified Bessel functions [13] with a = (d — 2)/2. From equation (f2"3~j) 
one obtains the asymptotic behaviour right at the transition point: 

N s ~ iV* , Co = Co* , N - oo , (24) 

with 

Co " { -(d-2) , d>2 • 
The crossover-exponent $ in equation ()24j) is given by JUj 

* = , 1 < d < 4 , d ^ 2 , (26) 

iV 5 ~ IniV for d = 2, N s ~ N/lnN for d = 4, and for d > 4 one finds $ = 1 
corresponding to a first-order transition. For Co < Co the finite fraction of adsorbed 
monomers scales like 

lim ^ ~ (Co - Co) K , Co < Co , (27) 

N— >oo iv 

with the exponent 

2- Id- 21 

K = | d _ 2 | ' 1 " d < 4 ' ^ 2 ' (28) 
For d = 2 one finds 

lim ^ ~ e~ 2 ^ , Co < , d = 2 , (29) 

N— ^oo iv 

which formally corresponds toK = oo. For d = 4 one finds lim N$/N ~ —1/ In (Co — Co) 
while for d > 4 the fraction tends to a finite value, which reflects the fact that in this case 
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the transition is of first order. Equations (J24|) - (|29|) are to be compared with equations 
© - (fT2"|) in section 11.11 where the variable Q — (q here corresponds to £ — £* there. 
In particular, for given (q < Co the cnam partition function x{to) exhibits a singularity 
for t \ ts(Co) where the function ts(Co) is determined by an analysis of the zero of the 
denominator in equation (|2"3*|) 10 j. The values of the exponents $ and k in equations 
f!26|) and ()2"%)) obey the scaling relation k = (1 — $)/$ from equation 1)12)1 . Note that for 
d > 2 the limit i? — > yields merely the trivial bulk result x( r ; ^o) — lAo> an d hence no 
phase transition for d > 2. Thus, in the present treatment it is necessary to keep the 
cylinder radius R finite even though the results for $ and k do not depend on R. 

Finally, we note that the adsorption of a Gaussian chain onto a rigid rod is 
equivalent to the denaturation transition of two flexible Gaussian chains A and B if 
monomer s of chain A can only interact with monomer s of chain B. This corresponds 
to an interaction of the form ~ J Q L ° ds 8\Ra(s) — Rb(s)] in the partition function (JT3J). It 
is easy to see that the system of two flexible Gaussian chains with the above interaction 
can be mapped on the system of one flexible Gaussian chain interacting with a rigid 
rod, using the transformation fTS) IT?) 

R(s) = R A (s) - R B (s) , Rcm(s) = ~ [R A (s) + R B (s)] ■ (30) 

For Gaussian chains (and only for Gaussian chains) the degrees of freedom described by 
the centre of mass (CM) coordinates Rcm(s) and the relative coordinates R(s) decouple 
from each other. The relative coordinates R(s) describe a Gaussian chain which interacts 
with the origin at R = in the d- dimensional subspace, while the degrees of freedom 
of the remaining D — d dimensions are unbounded and independent from the degrees 
of freedom of the d-dimensional subspace. By virtue of this mapping, the above results 
tflty - flZH) have also been obtained in [TH] . 

2.3. Renormalization of the field-theory bounded by the curved cylinder surface 

We now turn to the renormalization of the n- vector model defined by equation (|22[). 
The objective is to determine the scaling behaviour of the renormalized chain partition 
function Xren{R, t, c) in the grand canonical ensemble in terms of the cylinder radius R 
and renormalized parameters t (conjugate to the renormalized chain length L) and c. 
The n-vector model in equation can be dimensionally regularized and renormalized 
by minimal subtraction of poles in e = 4 — D. The renormalizations of the bulk field 
$(r), r G V, and the bulk parameters to, «o have the same form as in the unbounded 
case, and are given by [311133 E2] (we follow the convention of p?2") ) 

$ =zl /2 $ ren (31) 
t =fi 2 Z t t + t b (32) 

uo = h £ 2 d tv d/2 Z u u . (33) 

The parameters $ re n, t, u are renormalized counterparts of the bare parameters <3>, t , u 
in equation ((221) • I n a regularization scheme using a large momentum-cutoff A, the bulk 
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renormalization factors Z&, Z t , Z u absorb divergencies logarithmic in A, corresponding 
to poles in e in dimensional regularization. The parameter i& absorbs divergencies 
quadratic in A and describes the shift of the critical temperature T c of the n-vector 
model due to fluctuations. In dimensional regularization, i 6 = 0, and jSHl HUH E2| 



1 



n + 2 
36 e 



u + 0(u 3 



z f z 



t 



1 n + 2 

l + ^—u + 

3e 

3e 



(n + 2)(n + 5) n + 2 
9^2 6~7~ 

(n + 8) 2 



u 



+ 0(« 3 



9£ S 



3n + 14 

6e 



It 



0(w 3 



(34) 
(35) 

(36) 



The presence of the cylinder surface S requires, in addition, renormalization of the 

For a 

(37) 
(38) 



surface field $|s = $(r), r £ S, and of the surface parameter cq in equation 
flat surface these additional renormalizations are given by j3U E2J E2j 



$1 



7 x / 2 cb I 

£\ ^»ren|5 



Co = /i2 c c + c sp 

defining the renormalized surface field (<&\s)ren and surface parameter c. The new 
renormalization factors Z\ and Z c absorb divergencies logarithmic in A which occur at 
a flat surface, corresponding to poles in s in dimensional regularization. The parameter 
c sp absorbs divergencies linear in A and describes the shift of the multicritical point due 
to fluctuations (compare figure El — c ~ £). In dimensional regularization, c sp = 0, and 
Z x and Z c are given by jHHSl E2 

~(n + 2)(n + 5) n + 2' 



Zi 



n + 2 
l + ^—u + 

3e 

, n + 2 
3e 



9e 2 

(n + 2)(n + 5) 



+ 



35 

n + 2 



u 



+ 0(« 3 ) 



(39) 



4tt z 



u 2 + 0(u 3 ) (40) 



9e 2 ' 36e 

Equations ()37|) and (|3*H|) hold for a flat surface. As shown by McAvity and Osborn 
the required renormalizations become modified if the surface S is curved, like in the 
present case. While the renormalization of the surface field ($|s) ren remains unchanged, 
the renormalization of the surface parameter Co requires an additional, additive term 
that depends on the mean curvature [471] : 

liZ c c+ ^--^C(u,e) (41) 



or, with £ 



c 

Rcq and £ 
Co 



R 

fiRc, 

Z c ( + (d-l)C(u,e) . (42) 

Z c is the same renormalization factor as in equation ()38|) for a flat surface and C(u,e) 
to second order in u can be deduced from reference jlOj: 

n + 2 



C{u,e) 



9 



u 

- + 

e 



n + 5 n+1- 47r 2 



3 



12 £ 



u 



+ 0(u 6 



(43) 



Polymer adsorption transition 



13 



To proceed, we consider the two-point correlation function 
($ 1 (r s )$ 1 (r , )) = G(r s ,r / , i?;i ,Co,«o) , 



(44) 



where is located on the cylinder surface S and r' e V. From G, the chain partition 
function x(R> ^o, Co> u o) follows by an integration over r', see equation (J3TJ); by symmetry, 
X does not depend on rg. The renormalization group (RG) equation for the renormalized 
counterpart G ren of G follows in the standard way, using the relation 



G ren^S-, r> \ R] t, (, U 'i f) ~ 1 %1 ^ ^ G(r S , r' , R] t Q , ( , % 



(45) 



and the fact that G does not depend on //: fid^G] 



fi=0 



1 d( 
where we have used the abbreviation 



d f 



G , 



A^uU=o + P(u)d u - (2 + r) t )td t 



with 



(5{u) = fiB^ul^o 



0. This leads to the RG equation 
(46) 

(47) 
(48) 



and the exponent functions r)i(u) = fid^ In Zj| M=0 with i = $,1,1 The new feature of 
the RG equation generated by the surface curvature is the function 



jtt=0 



C + 



d — 1 n(d — 1) 



where 



Vc(u) 



/i(9 M lnZ c | M=0 



(49) 



(50) 



„ + V±l ( 4?r 2 _ l)u 2 + 0(u 3 } 



3 1* 

is the exponent function for c corresponding to a flat surface, with Z c from equation 
(|4()jl [3*2*] . A necessary condition for the two-point correlation function G to be scale 
invariant (SI) is that the right hand side of equation (j4Uj) vanishes, which is either the 
case for u = (Gaussian model) or if ( takes the value for which the square bracket in 
equation vanishes: 

C = Csi = -^ + 0( M ). (51) 

Note that the value of £si to leading order in u, Csi = ~(d — l)/3, is different from 
the value Q in equation (|23j) corresponding to the onset of the adsorption transition 
in the Gaussian model. It is also interesting to note that for a sphere, corresponding 
to d = D = 4 — e, the value of ("si to leading order in u, £si = —1, coincides with the 
value (ci — —(d — 2)/2 for which the Gaussian two-point correlation function at t = 
is conformal invariant (CI) |4*6l H7j. However, within the <3> 4 -model the special value £si 
is already fixed by scale invariance, whereas in the Gaussian model the value £ci is only 
fixed if one requires conformal invariance. 
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By solving the RG equation (J46j) for the two-point function in the standard way 
using characteristics and integrating over r', one derives the general scaling behaviour 
of the chain partition function in the grand canonical ensemble: 

Xren(R, t, C, «! ~ © (a^Rt v , (3^-^) (52) 

where 

AC = C - Csi (53) 
with Csi from equation ()51|). The exponent v — (2 + r^* ) _1 is a bulk critical exponent, 
while $fl at = v(l + 77*) (the crossover exponent) and 71 = u(2 — rj$ — rjl/2) are 
critical exponents associated with a flat surface [32] ■ The exponent functions rji(u) 
for i = t,c, $, 1 are defined below equation (|4*%|) and in equation (jBTlj) . The values 77* 
are the values of the exponent functions at the fixed point u*. The constants a and /3 
in equation ()52)1 are nonuniversal prefactors, while the function Q(x,y) is a universal 
scaling function. 

Equation (p^j) is in line with the Gaussian model (where v = $ flat = 1/2 and 
7i = 1): 

1 1 /2 

x(R, to, Co) ~ ©o(-R% j Co) j Gaussian model , (54) 
compare equation ((221) • Finally, equation (p)2"|) can be compared with the corresponding 
behaviour for a flat surface: 

Xren(t, c, m; /i) ~ t~ 71 ©fiat (/?'c £ _<1>flat ) , flat surface . (55) 

Note that this limit can be obtained from equation ()52)1 by rewriting the second scaling 
argument as AC^ - *^ = Rt u ■ Act~* flat . Equation (|B3|) then follows as the limit R — > 00 
of the scaling form Xren ~ ©(«' >-Rt" ', P' 'Act _ * flat ), where Ac = c + 0(l/i2). 

Let us come back to equation (|52*|). According to recent estimates one has 
$ flat = 0.518 jH| and v = 0.588 (HSJ for n = in £> = 3 so that the exponent v - $ flat 
in equation ()52|) is small but positive. From a naive point of view this would imply that 
the scaling variable AC in the second scaling argument of is irrelevant and could be 
omitted from the outset; however, one should keep in mind that the radius R in the first 
scaling argument is also irrelevant in principle. Now, the relevant question is whether the 
scaling function 6(2, y) exhibits a singularity on a certain subset of (x, y), corresponding 
to the polymer adsorption transition; compare the related discussion below equation © 
and below equation (J29|) . In fact, we expect this singularity to occur for y = ys(x) < 0, 
where now ys (x) is a function of the first scaling variable x = a[iRt v . Note that the 
scaling function 6 (x, y) in equation (JSIJ), corresponding to the Gaussian model, exhibits 
this kind of singularity indeed; compare section 12.21 Thus, in the present description, 
the adsorption transition is characterized by a balance of the two irrelevant variables 
AC and R. In this sense the scaling variable AC can be considered as a dangerously 
irrelevant variable. 

Finally we note that equation (p)2*|) implies a corresponding scaling form of the chain 
partition function with fixed chain length L: 

Z ren (R, L, C, u- n) ~ L 11 ' 1 ^ (afiRL- v , ^AC^~* flat) ) , (56) 
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with (different) nonuniversal prefactors a and (3 and a universal scaling function 
From equations (|52|) and ()56|) the number of adsorbed monomers Ng for A£ = and 
the finite fraction of adsorbed monomers N$/N for A£ < can be derived as outlined 
in section [TTTl in principle. However, the thin rod limit R — > corresponds to singular 
limits of the scaling functions and \1/ which are rather difficult to obtain. At least it is 
easy to see that the exponents $ and k defined in equations © and (fTTj) are universal, 
using the fact that and \I/ are universal scaling functions. To proceed, in the next 
section we use a different method to obtain estimates for the exponents $ and k for a 
cylinder in D = 3. 

2. 4- Estimates for the exponents $ and k by using the additivity of co-dimensions 

In this section we obtain estimates for the exponents $ and k for a cylinder in D = 3 
introduced in equations © and (fTTj) by means of an interpolation procedure between 
two known cases. Firstly, in the Gaussian model one has $ = 1/2 for d — 3, see 
equation corresponding to the point (d,D) = (3,4) in figure EJ Likewise, k = 1 
from the scaling relation k = (1 — $)/$, or from equation (|28|). for (d,D) = (3,4). 
Secondly, one has $ = on the whole line v~ 1 (D) — d = in figure 0] This result can be 
obtained by using the co-dimension additivity theorem, stating that the co-dimension of 
the intersection points of two objects of dimensions D\ and D 2 is given by the sum of 
their co- dimensions: D — D int — (D — D x ) + (D — D 2 ), i.e., 

D mt = D 1 + D 2 -D . (57) 

For example, two-dimensional surfaces generically intersect along curves in D = 3 
(Dint = 2 + 2 — 3 = 1) and only at isolated points in D = 4 (D int = 0). Equation 
(157}) is also expected to hold if one or both objects are fractal. In the present case, one 
object is a self-avoiding random walk with fractal (Hausdorff) dimension v ~ 1 and the 
other one is a "generalized cylinder" with co-dimension d (see figure EJ); the dimension 
of intersection points of these two objects is thus given by 

D int = v-\D)-d . (58) 

In figure IH the line Di nt = as a function of d and D is shown as the blue dashed line. 
An unbounded, free, self-avoiding random walk does not intersect with "generalized 
cylinders" located above the dashed line, apart from exceptional cases. In this sense, 
"generalized cylinders" above the dashed line are irrelevant perturbations for a free, self- 
avoiding random walk. Now, "generalized cylinders" located right on the line D int = 
correspond to marginal cases: An unbounded, free, self-avoiding random walk does 
intersect with "generalized cylinders" located on the dashed line, but only at isolated 
points. We thus expect that the number of intersecting monomers N$ grows with N 
for iV — > 00, but only logarithmically, i.e., N$ ~ In AT, which implies $ = 0; compare 
the case d = 2 for the Gaussian model discussed in section I2.2| and compare the case 
£ = £* with $ = in section 11.11 It should be noted that this argument only applies 
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to unperturbed random walks, and does not make any statement for walks that interact 
with the body. 

Thus, the values of the exponent $ at the end points of the green line in figure |U 
are available. This can be used to obtain an estimate for $ for an ordinary cylinder in 
D = 3 as follows. The shape of the dashed line in figure |U is known quite accurately 
by means of the e-expansion of v(D) in conjunction with the exact value v = 3/4 for 
n = in D = 2 |3*Hl 13*9"] . Thus, one may estimate $ for a cylinder in D = 3, located 
at the point (d, D) = (2, 3) in figure EJ by means of a linear interpolation between the 
known values of $ at the end points of the green line (compare references In 
this way we find for a cylinder in D = 3 the estimates, using equation (fT2"|) . 

1 1 - $ 

$ ~ - , k = — - — ~ 5 . (59) 
6 Q 

Since these exponents are universal and do not depend on the cylinder radius R, they 
are also expected to hold for a rigid rod with vanishing radius, or for a line of lattice 
sites in a numerical simulation of this system. 



3. Conclusion 



We have investigated the adsorption transition of a long flexible self-avoiding polymer 
chain onto a rigid thin rod by field-theoretical methods. The rod is endowed with a 
short-ranged adsorption energy S for the chain monomers so that, on increasing S, at 
some threshold value S* the chain undergoes a transition from an unbound state to 
a bound state, as shown in figure ^ The main results and remaining questions are 
summarized below. 

1) By means of general scaling arguments we obtained the scaling relation ()12)) for the 
exponents $ and k defined in equations (j3J) and (JUJ), and the phase diagrams shown 
in figure El in terms of the number of chain monomers N, the chemical potential fi 
conjugate to N, and the adsorption energy £. 

2) By representing the rod by a cylinder of finite radius R we could use available results 
for field theories with curved boundaries jlU]; see figure El By using renormalization 
group arguments, we derived the scaling behaviour of the chain partition function 
in the grand canonical ensemble, equation (j52j) . and in the canonical ensemble, 
equation (J56*l) . where L ~ N and t ~ /i. Notable features of the scaling results are 
the distinct form of the scaling variable ( oc Rc, where the parameter c is related 
to the surface potential for chain monomers, and the curvature-induced shift of £ 
in equation with £si from equation (|5*T|) . It also follows that the exponents $ 
and n introduced in equations ® and (fTTj) are universal. 

3) Because the cylinder radius R enters the scaling functions and \1/ in equations ()52j) 
and (foT)|) explicitly it is difficult to obtain the universal exponents $ and k directly 
from them. Therefore we used the co-dimension additivity theorem in conjunction 
with an interpolation procedure, as shown in figure |3J to obtain the estimates for 
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$ and k in equation (|59|). The check of these exponents and the scaling relation 
()12j) is a possible starting point for numerical simulations of this system. 

4) It would be interesting to introduce new methods to derive the exponents $ and k,, 
possibly avoiding the introduction of a finite cylinder radius R from the outset. 

5) It would also be interesting to explain the relation between Q and Csi discussed 
below equation (|5Tj) . 
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Appendix A. Mapping of the polymer system on the n-vector model 



The p 2 (r) interaction in equation (|T3~j) can be linearized by means of a Gaussian 
transformation [SZ1 E] • This procedure makes use of the Gaussian integral 



VX exp 



--X 1 AX + V X 
2 



, A 
det — 
2tt 



-1/2 



exp 



-b T A-H 

2 



(A.1) 



where X is a vector with discrete or continuous indices and the symmetric matrix A 

must have a positive definite real part. Using X(r) oc io~(r) with purely imaginary cr(r), 

3 

the matrix A(r, r') = — S(r — r'), and b(r) = ip(r), one finds 

d D rp 2 (v) 



exp 



3 
u 
u 
6 



(A.2) 



J Do" exp ^- J d D ra 2 (r) - J d D ra(v)p(v) 



Note that A is positive definite due to our assumption u > 0. Inserting (jA.2|) in (|T3^) 
yields 



Z (2) (r,r';L) = / Va exp |- ^ d D ra 2 (r) 



(A.3) 



x 



PR exp 



L dS (^) ~ ' d ° r [1/(r) + a(r)] P(r) 



The p 2 (r) interaction in equation (|13p has been replaced by the interaction of p(r) 
with an external, fluctuating potential er(r). The second line of equation (|A.3|) can be 
interpreted as the path integral representation of the evolution operator (r' '\e~ LH \v) 
in imaginary time s of a quantum-mechanical particle with Hamiltonian H = — A + 
V(r) + cr(r). The Laplace transform of this evolution operator with respect to L yields 
the resolvent 

1 ^ (A.J) 



/>oo 

/ dLe- tL (r' 


e- Lk 


r) = /r' 


Jo 







-A + 1 + V(r) + cr(r) 
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The resolvent can be represented in the standard way by the two-point correlation 
function of an n-component field <3> = ($1, . . . , $„) in the limit n — > 0. The result is 

" 3 

2u 



r(2)/ 



r,r'; L) 



Pcrexp 



(A.5) 



x C t ^L lim / P$$i(r)$i(r')e" s{ * } 

n->0 J 



with the action 

S{$} 



d D r 



(A.6) 



The Gaussian integration in equation (jA.5|) can be carried out using equation (|A.1|) 

3 

with the same X(r) oc icr(r) and matrix A(r, r') = —5(r — r') as before, and now 

u 



b{r) 



$ 2 (r). This leads to equations flEJ) - (fTTIl . 
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